Searching for dark matter isocurvature initial conditions with N-body Simulations 
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Small fraction of isocurvature perturbations may exist and correlate with adiabatic perturbations 
in the primordial perturbations. Naively switching off isocurvature perturbations may lead to biased 
results. We study the effect of dark matter isocurvature on the structure formation through N-body 
simulations. From the best fit values, we run four sets of simulation with different initial conditions 
and different box sizes. We find that, if the fraction of dark matter isocurvature is small, we can not 
detect its signal through matter power spectrum and two point correlation function with large scale 
survey. However, the halo mass function can give an obvious signal. Compared to 5% difference on 
matter power spectrum, it can get 37% at 2 = 3 on halo mass function. This indicates that future 
high precise cluster count experiment can give stringent constraints on dark matter isocurvature 
perturbations. 



I. INTRODUCTION 



Recent high accuracy observations, such as cosmic mi- 
crowave background radiation (CMB) 1], and the large 
scale structure (LSS)[2|, provide us wealthy information 
on the universe. Currently, so called concordance cos- 
mological model, in which approximately scale-invariant, 
Gaussian, adiabatic primordial perturbations seed the 
structure of our universe, is mostly favored. However, 
there exist some models which predict non-negligible 
isocurvature perturbations, such as multiple scalar fields 
inflation [3l-(l3| and curvaton |14| - ll6| models. The primor- 
dial isocurvature perturbations then can induce the cos- 
mological matter perturbations such as cold dark matter 
(CDM), baryons, dark energy and neutrino [l 7.-.20.] in the 
radiation era. 

Although the pure isocurvature perturbations have al- 
ready been ruled out by the observation of Boomerang 
and MAXIMA- 1 [2]| . models with primordial perturba- 
tions comprised of dominate adiabatic and a small frac- 
tion of isocurvature modes |22h34J still survive. Due to 
the quality of data, we still can not confirm the destiny of 
the isocurvature perturbation. Moreover, if we switch off 
the isocurvature perturbations naively, the result would 
be misleading (201 ISa] . So, it is important to seek a way to 
detect the isocurvature perturbations and give stringent 
constraints on isocurvature perturbations parameters. 

In Ref. [3J|, we have given constraints on parameters 
of two different cosmological models with latest astro- 
nomical data. One is with pure adiabatic initial condi- 
tion (IC) and the other is with mixed IC in which small 
fraction of isocurvature mode is correlated with dominat- 
ing adiabatic mode through cos A. We found that com- 
pared to adiabatic mode, the isocurvature perturbations 
have smaller amplitude {Af°/A'^'''^ ~ 10~^) and bluer 
tilt (n**° ^ 2). That is to say the isocurvature modes is 
pronounced on small scales (/c > 1 h Mpc"""^). However, 
on these scales, the structure grows non-linear ly, so we 
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can not study the details of evolution in a full analyti- 
cal way and then test them against future high-precision 
observations. 

On the other way, N-body simulations plays more and 
more important role in cosmology nowadays, especially 
in the study of nonlinearity. It can make predictions and 
compare with observations for a specific model. It can 
also be used to check the validity of a particular method 

In this work, we study the effects of CDM isocurvature 
perturbations on LSS by implementing N-body simula- 
tions. Given the best fit parameters in Ref. [3J], we 
carry out four sets of simulations with different ICs in 
different boxes. We find that, unlike the power spectrum 
and two point correlation function, mass function is sen- 
sitive to dark matter isocurvature perturbations. This 
indicates that, we can give stringent constraints on dark 
matter isocurvature perturbations with future high pre- 
cision cluster counts experiment. 

The structure of this paper is organized as follows. In 
Sec. ini we briefly introduce the CDM isocurvature per- 
turbation and set the ICs for N-body simulation. In Sec. 
Illli we describe the details of our N-body simulation and 
present the result. We give summary in Sec. HV] 



II. CORRELATED ADIABATIC AND 
ISOCURVATURE PERTURBATIONS 

Generally, to calculate matter power spectra -P(fc) nu- 
merically, one start from the time tic deep in the radiation 
dominated era, when all interesting scales of perturba- 
tions are outside the horizon. However, tic is different 
from the time t* when the corresponding mode k exits 
horizon during inflation. Therefore, One need transfer 
function Tij to transform perturbations from t* to tic- 

As we know, in the absence of isocurvature perturba- 
tions, the adiabatic perturbations TZ are conserved on 
superhorizon scales. On the contrary, the isocurvature 
perturbations S can evolve on superhorizon scale and 
can also seed adiabatic perturbations. Thus, the transfer 



function Tij can be written as|37| 
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where T^ is model dependent. To investigate isocur- 
vature perturbation without making use of any specific 
model, one often parametrizes the transfer function in 
the simple form of power law. 

For the Gaussian statistics, which is predicted by infla- 
tion, the power spectra characterize all the information. 
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where ^"1 = 7?. and X2 = S. We can parametrize primor- 
dial power spectra 'P{k) at tic as 
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where fcg is pivot scale, and both A^J and n'-J are 2 di- 
mensional symmetric matrices which denote the ampli- 
tude and power index, respectively. The amplitude of the 
spectra A^J can be written as 
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where cos A = ^adi,iso^^^adi^iso describes the correla- 
tion between adiabatic mode and isocurvature mode [6|, 
Af^^ and Af° stand for the amplitude of adiabatic and 
isocurvature modes, respectively. The spectra index nj' 
is 



^u - 
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with rij^' and nf° being the spectra indices for adiabatic 
and isocurvaure modes. Here, we have used the approx- 
imation n™"^ = 



'-i-< 



-i20] for simplicity. 



Since both adiabatic and isocurvature perturbations 
seed the large scale structure as 
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Then besides the normal adiabatic term, two other terms, 
isocurvature and cross-correlation terms, emerge in the 
expression of matter power spectrum, 

Af'P'"^\k) + A,'°P''°{k) 
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where P'(fc) can be described as 
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with T^{k) being transfer function of matter for IC i. 

We have given constraints on these parameters with 
latest observations in Ref. [^. The the best fit values 
are listed in Table HI Using CAMBJlil, We also sketch 
the power spectra in Fig. [T]with best fit values. 



TABLE I: The best fit value for models with different initial 
conditions. 



Parameters 


Aidabatic 


Mixed 


^6 


0.046 


0.044 


i Ijjl 


0.280 


0.267 


nA 


0.720 


0.733 


h 


0.700 


0.714 


lO^'Af' 


2.176 


2.420 


nf' 


0.960 


0.965 


lO^^Al"" 


- 


0.081 


n'f 


- 


2.716 


cos A 


- 


0.173 


0-8 


0.820 


0.865 




k(hMpc"') 



FIG. 1: Linear matter power spectra for the models in Table 
HI The red solid line corresponds to standard ACDA'I model 
with adiabatic initial condition while the blue dash-dotted 
line is given by mixed initial condition. 



III. N-BODY SIMULATION 

We perform our simulations with GADGET-2i[3il, 
a massively parallel TreePM-SPH (Tree Particle Mesh- 
Smoothed Particle Hydrodynamics) code. For collision- 
less particles, the gravitational field is calculated with 
a low-resolution particle-mesh(PM) algorithm on large 
scales, while forces are delivered by tree on small scales. 
We do not use the SPH part since only cold dark matter 
particles are considered in this work. 



A. Initial Conditions 

With power spectrum plotted in Fig. [l] we can gen- 
erate positions and velocity ICs for particles at cosmic 
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TABLE II: Simulation details for adiabatic and mixed initial 
conditions. 





Set I 


Set II 


Initial Condition 


adiabatic 


Mixed 


adiabatic 


Mixed 


Ltox{h~^Mpc) 

Npart 

Lsoft{h~'^kpc) 

Zstart 


1000 
512=* 

10 

49 


1000 
512^ 

10 

49 


100 

320^ 

5 

49 


100 

320^ 

5 

49 



time T using the Zel'dovich Approxiinations(ZA). 

x(q,T) = q + i?+(r)*(q), (9) 



v(q,T) = D+{T)^{q), 



(10) 



where x is the perturbed comoving coordinates and v = 
^ is the proper pecuhar velocity, q, the lagrangian co- 
ordinates generated from glass configuration [401, denote 
the unperturbed comoving position. D~^{t) is the linear 
growth factor normalized to z = 0. '4'(q) is displacement 
field calculated from the density fluctuation field which is 
the convolution of a random white noise with the square 
root of the linear power spectrum[41|. 

The ICs are set a,t z = 49 when the second order La- 
grangian perturbations correction can be ignored safely, 
we run four sets of simulations with different box sizes 
to explore the differences between two initial conditions 
on different scales. The larger boxes whose length is 
1000 h^^ Mpc provide good statistics on large scales from 
k ~ 10~^ h Mpc^^ to fc ■^ 1 h Mpc~^, while the smaller 
boxes with L = 100 h^ Mpc can give high resolution 
extending to fc ^^ 10 h Mpc" . In Set I, the mass reso- 
lution is 5.8 X 10" h"^Mo with Np = 512^ while in Set 
II the mass resolution is about 2.4 x 10^ h~^MQ. The 
force resolution is taken as ~ 0.5% of the mean particle 
interval (Tab. HH. 

Because we can get only one value of as for a specific 
survey and to cease the effect of different choices of cos- 
mological parameters, we renormalize the power spectra 
in Fig. [T]to the same trg = 0-8. This setting may make 
results present below not so obvious, however, what we 
are interested in is the relative difference, which is inde- 
pendent on the renormalization. 



We calculate the 2-point correlation function for Set 
I at z = with the pair-count estimator proposed by 
Landy & Szalay[43]: 
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where DD and RR are the autocorrelation function of the 
simulation particles and randomly sampled points respec- 
tively, DR is the cross-correlation between the data and 
random points. From Fig. [51 we can find that, the posi- 
tion and width of BAO from which H{z) and Da{Z) are 
extracted, are almost the same for the two different simu- 
lations. This result is reasonable in two aspects. Firstly, 
the BAO observation mainly depends on the background 
parameter, while has little to do with the origin of the 
perturbations; secondly, we set ICs with best fit param- 
eters. The behavior on large scales {k < 0.2 h Mpc~^) is 
well constrained |44|. especially the position of first peak 
in the CMB angular power spectrum[l|. It implies that, 
if the isocurvature fraction is small enough, we can not 
discriminate two initial conditions from BAO observa- 
tion. This is an important systematic error in constrain- 
ing dark energy from BAO data|33l[35|. 
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FIG. 2: The 2-point correlation function for two simulations 
in the large box at z = 0. The blue stars denote mixed case 
while the red points are for adiabatic IC. 



B. Numerical Results 

1. Correlation Function 

The baryon acoustic oscillation (BAO), as a standard 
ruler, is a powerful tool to study the dark energy. It is 
also a useful tool to detect the dark matter isocurvature 
perturbation. The presence of dark matter isocurvature 
perturbation would alter the position of first peak in the 
CMB angular power spectrum, which is the right scale of 
BAO. The peak in the two point correlation function and 
the wiggles in the power spectrum are the useful tools to 
track the behavior of BAOlij. 



2. Power Spectra 

Power spectrum, defined as the Fourier transforma- 
tion of two point correlation function, is the key physical 
quantity in understanding clustering properties. With a 
Gaussian initial condition as selected in this paper, power 
spectrum gives a complete statistical description of fluc- 
tuations. 

'POWMES' is a power spectrum estimator based on 
the Taylor expansion of the trigonometric functions j45|. 
The further 'foldings' scheme makes it possible to give an 
accuracy measurement of power spectrum up to a scale 
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is the Nyquist frequency and n 



fold 



where fc„„ — -j ^ 

is the number of foldings which is set as 2 in this work. 
We plot the power spectra at different redshifts as well 
as their ratios of Set I in Fig. [3l 




FIG. 3: Top panehThe power spectrum of simulation Set I. 
From bottom to top: z = 5, z = 3, z = 1, z = 0; The blue 
stars stand for mixed initial condition while red points denote 
adiabatic ACDM. Bottom panel: the ratio of matter power 
spectrum between two initial conditions. The thickness is 
proportional to the scale factor. 

From Fig. [S] we can find that power spectra of sim- 
ulation with mixed IC is smaller than the one in adi- 
abatic case on large scales (small k). r, the ratio of 
Pmixed to Padiabatic, grows as & function of time on large 
scales, A: > 2 h Mpc^^ and decreases on small scales, 
A: < 2 h Mpc"^ For k ~ 0.02 h Mpc"\ r reaches about 
0.95 at ^ = 0. That is to say, the largest discrepancy in 
power spectrum is about 5%. 

What we have to keep in mind is that we have 
renormalized the initial power spectra to the same erg 
and the weight fc^ C^'^^J — ^^^jpp-Y peaks around k ~ 
0.3 h Mpc^^. So, the power spectra we got from N-body 
simulations with mixed IC should be similar to the one 
with adiabatic IC on large scales while be larger on small 
scales without renormalization. 



3. Halo Mass Function 

Mass function is defined as the abundance of dark mat- 
ter haloes in a specific mass ranges. It is a key quantity to 
describe the large scale structure in the nonlinear regime. 
Press and Schechter firstly provided the theoretical de- 
scription in a simple spherical collapse model[47|. Sub- 
sequently, people made some improvements on this sim- 
ple modelling and introduced more complex ellipsoidal 
collapse models [4a]. Meanwhile, a lot of literatures try 
to fit the halo mass function in the manner of N-body 
simulation EiJii. 



We identify haloes with AHF^ (Amiga's Halo 
Find) [46], an adaptive mesh based finder. After placing 
grid across the box and further refinement, AHF assigns 
each particle to a grid with cloud-in-cell (CIC) interpo- 
lation. Then, AHF probes the halo at each density peak 
using spherical overdensity (SO) algorithm. The radius 
of sphere is grown until the interior density reaches a 
specific value 
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where pbkg = ^mPcriti^ + z)^ is the mean density of 
whole box. To compare mass functions of these two dif- 
ferent models, we set A as 200 in this paper. With these 
scheme, AHF can find all structures and substructures 
simultaneously. Moreover, we only keep haloes with at 
least 20 dark matter particles, i.e. the mimimum mass 
of haloes is around 4 x 10^°/i~^Mq. 

We introduce the cumulative mass function which is 
defined as mean number densities of haloes with mass 
larger than a specific mass, 
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where N{> M) is the number of haloes with mass greater 
than M. This is related to the mass function f{(j) through 
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The cumulative mass function as well as the ratio r = 
nmixed/nadia for Simulation Set II at redshifts z = 3, 1, 
are sketched in Fig. S) We can find that the mass 
function is almost the same on large mass scale, from 
lO^^Af0 to 1O^^M0. However, there are more haloes for 
adiabatic IC than for mixed IC on small mass scale ( 
M < 1O^^M0 ). Moreover, the discrepancy increases 
as the mass gets smaller. For M ~ 5 x IO-'^^Mq, the 
difference is about 26% at z = 0. We also plot the ratio r 
against redshift z for M = 5 x 10^°/i~^AfQ in the bottom 
right of Fig. H) We find that the ratio decreases as time 
evolves. Compared to r = 74% at z = 0, the ratio at 
z = 3 is only 63%. This behavior can be ascribed to the 
late-time non-linear evolution. That is to say, a high- 
redshift survey is helpful to seek the DM isocurvature 
perturbations signal. 

Since we have renormalized to the same erg , the power 
on large scales for mixed IC is less than for adiabatic IC 
(Fig. [3]), and there should be more haloes for mixed IC 
than for the adiabatic one without renormalization. 
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FIG. 4: Cumulative mass function as well as the ratio of two situations for simulation Set II. Top left are for z — 3, top right 
z — 1 and bottom left z — 0. The bottom right is the plot of discrepancy evolvement at M = 5 x W^'^Mq. 



IV. SUMMARY 

Isocurvature perturbations, inevitably generated from 
multi-field inflation or curvaton models, can be used to 
test these models. Although pure isocurvature pertur- 
bation models have already been ruled out, there still 
exists possibility that a small fraction of isocurvature 
mode is correlated to the dominating adiabatic pertur- 
bation. This is important to the parameter estimation, 
since rough ignorance would lead to biased result. 

With the best fit values obtained in Ref.jSJ], we per- 
form four sets of N-body simulations to seek a best way 
to detect the isocurvature perturbations. We find that, 
if the fraction is small enough.{Aiso / Aadia ^ 3%) we can 
not peek it in the BAO observation. The position and 
the width of the bump in two-point correlation function, 
which mainly depend on the background parameter, are 
almost the same. There are some differences in matter 



power spectrum and halo mass function. However, the 
5% difference in matter power spectra makes it hard to 
be observed. On the contrary, the deviation in the halo 
mass function is obvious. The difference is getting larger 
as we go to higher redshift. For M ~ 5 x lO^^M© the dis- 
crepancy can get 37% at z = 3. This implies that, with 
future precise cluster number count observations, we can 
detect the initial condition of dark matter isocurvature 
and give stringent constraints. 



Acknowrledgements 

We perform our simulations on Deepcomp7000 of Su- 
percomputing Center, Computer Network Information 
Center of Chinese Academy of Sciences. We thank Shi 
Shao, Wenting Wang, Jia-Xin Han, Hong Li, Jun-Qing 
Xia, Yi-Fu Cai Taotao Qiu and Mingzhe Li for helpful 



discussions. 



[1] E. Komatsu et at [WMAP Collaboration], Astrophys. J. 

Suppl. 192, 18 (2011). [28 

[2] B. A. Reid et al, Mon. Not. Roy. Astron. Soc. 404, 60 

(2010). [29 

[3] A. R. Liddle, A. Mazumdar, F. E. Schunck, Phys. Rev. [30 

D58, 061301 (1998). 
[4] R Kanti, K. A. Olive, Phys. Rev. D60, 043502 (1999). [31 

[5] E. J. Copeland, A. Mazumdar, N. J. Nunes, Phys. Rev. 

D60, 083506 (1999). [32 

[6] D. Langlois, Phys. Rev. D59, 123512 (1999). 
[7] D. Wands, N. Bartolo, S. Matarrese and A. Riotto, Phys. [33 

Rev. D 66, 043520 (2002). 
[8] Y. S. Piao, R. G. Cai, X. m. Zhang and Y. Z. Zhang, [34 

Phys. Rev. D 66, 121301 (2002). 
[9] S. Dimopoulos, S. Kachru, J. McGreevy and [35 

J. G. Wacker, JCAP 0808, 003 (2008). 
[10] D. Langlois, S. Renaux-Petel, D. A. Steer and T. Tanaka, [36 

Phys. Rev. Lett. 101, 061301 (2008). [37 

[11] Y. -F. Cai, H. -Y. Xia, Phys. Lett. B677, 226-234 (2009). 
[12] Y. -F. Cai, W. Xue, Phys. Lett. B680, 395-398 (2009). [38 

[13] Y. -F. Cai, J. B. Dent, D. A. Easson, [arXiv:1011.4074 [39 

[hep-th]]. 

[14] D. H. Lyth and D. Wands, Phys. Lett. B 524, 5 (2002). [40 

[15] T. Moroi and T. Takahashi, Phys. Lett. B 522, 215 [41 

(2001) [Erratum-ibid. B 539, 303 (2002)]. [42 

[16] D. H. Lyth, C. Ungarelh and D. Wands, Phys. Rev. D 

67 (2003) 023503. [43 

[17] D. Seckel and M. S. Turner, Phys. Rev. D 32, 3178 

(1985). [44' 

[18] A. D. Linde, Phys. Lett. B 259, 38 (1991). 
[19] A. D. Linde and V. F. Mukhanov, Phys. Rev. D 56, 535 [45 

(1997). 
[20] J. Liu, M. Li and X. Zhang, arXiv:1011.6146 [astro- [46 

ph.CO]. 
[21] K. Enqvist, H. Kurki-Suonio and J. Valiviita, Phys. Rev. [47 

D 62, 103003 (2000). 
[22] C. Gordon, A. Lewis, Phys. Rev. D67, 123513 (2003). [ 

[23] P. Crotty, J. Garcia-Bellido, J. Lesgourgues et al., Phys. 

Rev. Lett. 91, 171301 (2003). [49 

[24] M. Bucher, J. Dunkley, P. G. Ferreira et al, Phys. Rev. 

Lett. 93, 081301 (2004). [50 

[25] K. Moodley, M. Bucher, J. Dunkley et al, Phys. Rev. 

D70, 103520 (2004). [51 

[26] H. Kurki-Suonio, V. Muhonen, J. Valiviita, Phys. Rev. 

D71, 063005 (2005). [52 

[27] M. Beltran, J. Garcia-Bellido, J. Lesgourgues et al., Phys. 



Rev. D71, 063532 (2005). 

R. Bean, J. Dunkley, E. Pierpaoli, Phys. Rev. D74, 

063503 (2006). 

R. Trotta, Mon. Not. Roy. Astron. Soc. 375, L26 (2007). 

L SoUom, A. ChaUinor and M. P. Hobson, Phys. Rev. D 

79, 123521 (2009). 

J. Valiviita and T. Giannantonio, Phys. Rev. D 80, 

123516 (2009). 

M. Beltran, J. Garcia-Bellido, J. Lesgourgues and 

M. Viel, Phys. Rev. D 72, 103515 (2005). 

A. MangiUi, L. Verde and M. Beltran, JCAP 1010, 009 
(2010). 

H. Li, J. Liu, J. Q. Xia and Y. F. Cai, arXiv:1012.2511 
[astro-ph.CO]. 

C. Zunckel, P. Okouma, S. M. Kasanda, K. Moodley and 

B. A. Bassett, Phys. Lett. B 696, 433 (2011). 
J. S. Bagla, Curr. Sci. 88, 1088 (2005). 

L. Amendola, C. Gordon, D. Wands and M. Sasaki, Phys. 
Rev. Lett. 88, 211302 (2002). 

A. Lewis and S. Bridle, Phys. Rev. D 66, 103511 (2002). 
V. Springel, Mon. Not. Roy. Astron. Soc. 364, 1105 
(2005). 

S. D. M. White, arXiv:astro-ph/9410043. 
E. Sirko, Astrophys. J. 634, 728 (2005). 

D. J. Eisenstein et al. [SDSS Collaboration], Astrophys. 
J. 633, 560 (2005). 

S. D. Landy and A. S. Szalay, Astrophys. J. 412, 64 
(1993). 

B. A. Reid et al, Mon. Not. Roy. Astron. Soc. 404, 60 
(2010). 

S. Colombi, A. H. Jaffe, D. Novikov and C. Pichon, 

arXiv:0811.0313 [astro-ph]. 

S. R. KnoUmann and A. Knebe, Astrophys. J. Suppl. 

182, 608 (2009). 

W. H. Press and P. Schechter, Astrophys. J. 187 (1974) 

425. 

R. K. Sheth, H. J. Mo and G. Tormen, Mon. Not. Roy. 

Astron. Soc. 323, 1 (2001). 

R. K. Sheth and G. Tormen, Mon. Not. Roy. Astron. Soc. 

308, 119 (1999). 

A. Jenkins et al, Mon. Not. Roy. Astron. Soc. 321, 372 

(2001). 

M. S. Warren, K. Abazajian, D. E. Holz and L. Teodoro, 

Astrophys. J. 646, 881 (2006). 

J. L. Tinker et al, Astrophys. J. 688, 709 (2008). 



